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Abstract. We analyse a special case of the robust stabilization problem under struc- 
tured uncertainty. We obtain a new criterion for the solvability of the spectral Nevanlinna- 
Pick problem, which is a special case of the /^-synthesis problem of H°° control in which 
/i is the spectral radius. Given n distinct points Ai, . . . , X n in the unit disc and 2x2 
nonscalar complex matrices Wi, . . . , W„, the problem is to determine whether there is 
an analytic 2x2 matrix function F on the disc such that F(Xj) = Wj for each j and the 
supremum of the spectral radius of F(\) is less than 1 for A in the disc. The condition is 
that the minimum of a quadratic function of pairs of positive 3n-square matrices subject 
to certain linear matrix inequalities in the data be attained and be zero. 



1. Introduction 

In this paper we study an optimisation problem that arises in the design of a stabilizing 
controller for a linear time-invariant system that is subject to structured uncertainty. We 
show that, in a special case, the existence of a robustly stabilizing controller is equivalent 
to the condition that the minimum of a quadratic objective function of a matrix pair 
subject to a linear matrix inequality (LMI) be attained and be zero. 

Robust control theory provides a rigorous framework for the formulation and analysis 
of specifications of control systems for plants that are subject to sundry types of uncer- 
tainty; see for example [151 E] or [IS Chapter 8] . One of the tools of the theory is the 
structured singular value of an operator or matrix corresponding to a given uncertainty 
class ([15] or [18^ Definition 8.13]); this is a cost function that generalizes the operator 
norm. It is denoted by //, and leads to the "/i-synthesis problem", which is a problem 
of optimization over a class of analytic matrix functions in a disc or half-plane. Special 
cases of the //-synthesis problem are the Nehari and Nevanlinna-Pick problems, which have 
classical solutions, but in virtually no other case is there an analytic solution - there are 
only approximate numerical methods that have neither guaranteed convergence nor error 
bounds. 

The //-synthesis problem is an interpolation problem for analytic matrix functions. It 
is a familiar fact that robust stabilization leads to interpolation problems (|2CH 118]). For 
a given nominal plant and uncertainty class, the set of all stable closed-loop transfer 
functions can be parametrised, resulting in a class of analytic matrix functions F that 
are subject to interpolation conditions. To maximize the uncertainty region about the 
nominal plant that can be simultaneously stabilized one must minimize over F the quantity 
sup^ /j,(F(X)), where A varies over a disc or halfplane and fj, is a cost function that encodes 
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structural properties of the uncertainty set. In general the interpolation conditions that 
F satisfies are of the "model matching" type (Proposition 12.21 below), but in this paper 
we restrict ourself to the case that the values of F are specified at finitely many points. 

The computation of fi for an arbitrary block-structured uncertainty is known to be 
NP-hard |13( 127], However, one familiar and easily-computed instance of the structured 
singular value is the spectral radius r(-) of a matrix; it corresponds to a one-dimensional 
uncertainty class. Accordingly the following is an instance of the /i-synthesis problem. We 
denote the unit disc of the complex plane C by D. 

Problem SNP Given distinct points A\,...,\ n G D and target matrices W\, . . . , W n of 
type k x k find an analytic k x k-matrix-valued function F such that 

F(\j) = Wj for j = 1, . . . , n, and 

supr(F(A)) is minimized. 
AeB 

This problem is called the spectral Nevanlinna-Pick problem and has attracted attention 
as a test problem for //-synthesis ( [9j [TO] and other papers by these authors) . There are 
heuristic algorithms that calculate approximate solutions to //-synthesis problems ( \18\ 
Section 9.3.3] or |23|). These algorithms are in current industrial use, but they are nev- 
ertheless widely regarded as not fully satisfactory. In the absence of an adequate analytic 
theory it is even difficult to test their outputs for closeness to optimality. 

In this paper we prove a solvability criterion for the spectral Nevanlinna-Pick problem 
in the 2x2 case (that is, k = 2). The main result is the following. 

Theorem 1.1. Let \%, . . . , X n be distinct points in B and let W\, . . . , W n be 2x2 complex 
matrices, none of them a scalar multiple of the identity. Let Sj = tr Wj , Pj = det Wj for 
each j and let z±,Z2,z^ be any three distinct points in D. The following three conditions 
are equivalent. 

(1) There exists an analytic 2x2 matrix function F in B such that 

F(^) = Wj for j = l,...,n (1.1) 

and 

r(F(X)) < 1 for all A G B; (1.2) 

(2) the semidefinite program 

min (trAO 2 - tr(iV 2 ) (1.3) 



subject to 



N=[N l£ , ]k ]lf =1Ak=1 >0, 
M=[M iejk ]lf =iek=1 >0, 

> [(1 - z t z k )N U jk] + [(1 - AiAjJMtfjfc] (1.4) 



_ / 2zipi - s A 2z k pj - Sj 
V 2 - z e si J 2 - z k Sj 

is feasible and the minimum (jl.3p is attained and has value zero; 
(3) the semidefinite program in (2) admits a feasible pair (N, M) such that rank N < 1. 
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Here the inequality sign denotes the usual partial order on the space of 3ra-square 
Hermitian matrices. This result is a part of Theorems 18.11 and 18.41 which are established 
in Section [H 

Neither (2) nor (3) provides a convex program to resolve the existence of the desired 
function F in (1). In (2) the feasible region is convex, but the objective function (tr N) 2 — 
tr(iV 2 ), though non-negative and quadratic, is not concave, so that it may fail to attain a 
minimum at an extreme point of the feasible region and it may have many local minima. 
Condition (3) is a simple reformulation of (2), since it is easy to show that (triV) 2 — 
tr(iV 2 ) = if and only if rankiV < 1 (Lemma 18.21 below). In (3) there is no objective 
function, but the rank constraint means that the feasible set is not convex. Thus neither 
formulation is an LMI problem, as treated for example in [12]. We do not know whether 
(2) or (3) can be the basis of an efficient numerical procedure. There are, however, many 
papers in the literature on the numerical solution of optimization and feasibility problems 
of types that include (2) and (3), and it may be that some of the methods proposed will 
provide (in this special case) an effective alternative to algorithms in the current literature 
such as Q3H EE1 [23]. Rank-constrained LMIs have been studied in [251 E3 ED [22] , while 
there are many algorithms for the optimization of smooth functions over convex sets [23]. 
In the programs (2) and (3) the feasible pair (N, M) may be constrained to lie in a set 
having a known prior bound (Proposition II 1 . H below) . 

Of course relaxation of the rank constraint in (3) yields a necessary condition for the 
solvability of a 2 x 2 spectral Nevanlinna-Pick problem in the form of the feasibility of a 
true LMI: see Corollary W% 

The paper is organized as follows. Section [2] describes the robust stabilization problem, 
presents a concrete example and outlines the reduction of the problem to a model matching 
problem. Section [3] describes the symmetrized bidisc T and its magic functions. It also 
describes the reduction of a 2 x 2 spectral Nevanlinna-Pick problem to an interpolation 
problem in the space Hol(B, V) of analytic functions from ID) to T. Section 0] gives a worked 
example which is more general than the one in Section[2l This example, though illustrative, 
is limited to the case of systems with only two right-half-plane poles, and so motivates the 
need to develop an alternative approach. To this end Section [5] presents a duality between 
the space Hol(B, T) and a subset of the Schur class S2 of the bidisc. In Sections [6] and [8] we 
use Hilbert space models for functions in £2 to obtain necessary and sufficient conditions 
for solvability of the interpolation problem in the space Hol(B, T). Section [7] presents 
another approach to the realization of the relevant functions in S2 ■ In Section [9] we give 
some matricial formulations of the solvability criteria for the 2x2 spectral Nevanlinna- 
Pick problem. In Section [10] the procedure for constructing interpolating functions in 
Hol(B,r) developed in Section [8] is shown to be general: in principle it yields all possible 
interpolants. Section [11] contains some remarks about the numerical implementation of 
our procedure. 

The closed unit disc in C will be denoted by B _ and the unit circle by T. 

The complex conjugate transpose of a matrix A will be written A*. The symbol / will 
denote an identity operator or an identity matrix, according to context. For a matrix A 
and a non-zero scalar A we shall sometimes write A/X as a synonym for A _1 A The right 
half plane {s : Re s > 0} will be denoted by H, and RH°° will be the space of real-rational 
matricial functions that are analytic and bounded on H; the type of the matrices will be 
understood from the context. 
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2. Robust stabilization 

The theory of the robust stabilization of a plant subject to structured uncertainty is a 
particularly elegant chapter of H°° control. It was developed in the 1980s and 1990s (see 
|15[ [16] and many references in [IB]); it is well described in [181 Chapters 8 and 9]. In 
this section we sketch the reduction of the robust stabilization problem to an optimization 
problem of 'model matching' type. Consider the system S in Figure 1: 



G 



K 



G 



Figure 1 : £ Figure 2 : G 

£ is a feedback system with uncertainty; A,G and K denote finite-dimensional linear 
time-invariant systems, identified with their transfer functions, which are real rational 
matrix- valued functions of the frequency variable s. The nominal plant 



G = [Gi 



(2.1) 



is given and is proper. Uncertainty is modelled by the assumption that the 'true plant' 
G A is given by Figure 2 for some A which belongs to a prescribed 'uncertainty set' A 
but is otherwise unknown. We shall say that G is robustly stabilizable with respect to A if 
there exists a stable controller K such that the configuration in Figure 1 is stable for all 
A € A. Mathematically the requirements are that 

(1) the controller K belong to RH°°, 

(2) the system in Figure 1 be well posed for all A 6 A and 

(3) the system in Figure 1 be stable for all A S A. 

A system is said to be well posed if the transfer functions between different branches of the 
interconnection are well defined - in particular, if all the inverses occurring in the transfer 
functions exist in the ring of square rational matrix functions of appropriate type. See 
[T8[ page 282] for a fuller discussion of well-posedness. 

Here is an instance of the robust stabilization problem. For 7 S H let 6 7 denote the 
stable allpass function of degree 1 (or Blaschke factor) 



6 7 (s) 



7 



Let /, g G RH°° be defined by 



/(*) 



4s + 3 



3s + 1 



s + 7' 



fa) 



s e 



(s + 3)(s + 5) 
3(s + l) 2 ' 



It may be verified that 



fh + gb 3 = l. 



(2.2) 



(2.3) 



(2.4) 
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As in [TH1 page 262], define the uncertainty set A^o by 

Ai, = {51 : \S\ < 1}. 

Example 2.1. Let 



(2.5) 



G= [G 



ij 



3 . 



6163 


61 + 563 


1 


9 


1 





1063 + /61 


26^ + fei6 3 


/61/63 


1 





61/63 


1 




1 


fl/61 


1/61 








1 


//6 3 


1 





1/63 


1 





I/63 








61/63 





63/61 





I/61 


63/61 






(2-6) 



regarded as partitioned into 2x2 blocks. Does there exist a robustly stabilizing controller 
K for the system S of Figure 1 with respect to the uncertainty set Ai ? 

The answer to this question is given in Proposition 14.41 on page 1151 
For a general plant G, evidently some stabilizability assumption on G is a prerequisite 
for the existence of the K we are seeking. We shall assume that G is stabilizable, that is, 
there exists a controller K such that the lower loop of Figure 1, which is the system M(K) 
shown in Figure 3 below, is internally stable. With this assumption, by [181 Lemma 5.4], 
K stabilizes G if and only if it stabilizes G33. However, stabilization of G does not imply 
stabilization of G A for a general A. 
In the notation of [181 page 196], 



M(K) = [M, 



'J 



2 

i,j=l 



clef 



S(G,K) 

Gn G12 
G21 G22 



+ 



G13 

G23 



K(l - G33K)- 1 [G31 G 



32 



(2.7) 



G 



K 



M(K) 



Figure 3 : M(K) Figure 4 : L 

The system E of Figure 1 can also be written as in Figure 4, and so it has transfer 
function 



def 



S(M(K), A) = M 22 + M 2 iA(l - MiiA)~ l Mi2, AG A. 



(2.8) 



For a fixed K, Theorems 8.22 and 9.8 of [H] assert that S(M(K), A) exists and belongs 
to RH°° for all A G A if and only if K stabilizes G33 and 

(2.9) 



sup / u(Mn(s), A) < 1. 

s&W 



Here /u(-, A) denotes the structured singular value of a matrix relative to the uncertainty 
set A [181 Definition 8.13]. For the present paper it is enough to assume that A is such 
that A) = r(-), the spectral radius. 
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Now the set of all stabilizing controllers of G33 has an elegant parametrization due 
originally to Youla \18\ Theorems 5.13 and 5.14]. Let G33 have the doubly coprime fac- 
torization 



G33 = NM~ 



M~ l N 



(2.10) 



" X 


—Y 




M 


Y 


-N 


M _ 




N 


X 



over RH°° where N, M, X, Y, N, M, X and Y belong to RH°° and satisfy 

/ (2.11) 
Proposition 5.10]). 

(2.12) 



(every proper real rational plant admits such a factorization: see 
Then the general stabilizing controller of G33 is given by 

K = (Y — MQ)(X - NQ)- 1 = (X- QN)~\Y - QM) 



for some Q £ RH°° such that X(oo) — N(oo)Q(oo) is invertible. Moreover, for K in 
equation (I2.12|) . 

M{K) = fi - T2QT3 



(2.13) 



where 



J- 1 
f 3 



11 


G12 


21 




13 


M, 


23 





G 

Gi 



M[G 31 G 



+ 



G13 

G23 



YM[G 31 G 32 ] 



(2.14) 



■ VI 



To summarize: 

Proposition 2.2. Let G = [Gjj]^ be a stabilizable plant and let G 33 have the doubly 
coprime factorization (I2.10|) - (12. lip over RH°° . Let the zero matrix belong to A C RH°° . 
There exists a controller K £ RH°° such that the system E of Figure 1 is internally stable 
for all A £ A if and only if there exists Q £ RH°° such that X(oo) — N(oo)Q(oo) is 
nonsingular and 

(2.15) 



where 



sup n((T 1 -T 2 QT 3 )(s), A) < 1 

sen 



Ti = Gn + G 13 YMG 3l , T 2 = G l3 M, T 3 = MG 



31. (2.16) 

Moreover the general robustly stabilizing controller of £ for the uncertainty set A is given 
by equation (|2.12p for some Q £ RH°° such that X(co) — N(oo)Q(oo) is nonsingular and 
inequality ()2. 15j) holds. 

The problem of whether there exists a function Q with these properties is called a model 
matching problem. 

Proof. Suppose there does exist a controller K that stabilizes £ for all A £ A. Since 
0/ £ A, in particular K stabilizes the plant 



G° 



G22 G23 
G32 G33 
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Hence K stabilizes G33, and so by the Youla parametrization there exists Q G RH°° such 
that X(oo) — N(cx>)Q (00) is nonsingular and K satisfies equations (|2.12p . Furthermore K 
satisfies the inequality (|2.9|) . In view of equation (12.131) . 



M11 = [I 0] M(K) 



T x - T 2 QT 3 



where 



[/ 0]fi 



[1 o]f : 



2- 



It follows from equations (|2.14p that Ti , T 2 , T3 are given by equations (|2.16p . Thus necessity 
holds in Proposition 12.21 To prove sufficiency one simply reverses the steps. □ 



By the example on [181 page 257] fi(-, A^o) = r(-), the spectral radius. Thus the robust 
stabilization problem for the nominal plant G with uncertainty set A^o reduces to the 
following. 

Find Q G RH°° such that X(po) — N(oo)Q(oo) is nonsingular and 

supr((Ti-T 2 QT 3 )( S ))<l. (2.17) 

sen 

Now specialise further to the case that T2 and T3 are scalar matrix functions and the zeros 
of T2T3 in EI are simple, say s\, . . . , s n £ H. Then the set of functions T\ — T2QT3 as Q 
ranges through RH°° is just the set of functions in RH°° that agree with T\ at s\, . . . , s n , 
and the question becomes whether there exists a function i 7 G RH°° such that 

F(sj) = ri(sy) for j = 1, . . . ,n and sup r(F{s)) < 1, (2.18) 

sen 

which (after the application of a Cayley transform) is an instance of Problem SNP. In the 
case that Gu and G33 are 2x2 matrix functions Problem SNP can be analysed with the 
aid of the theory of the symmetrized bidisc, as presented in the next section. 

3. The symmetrized bidisc and its magic functions 
The open and closed symmetrized bidiscs are the subsets 

Q = {(z + w,zw) : \z\ < 1, \w\ < 1} (3.1) 

and 

T = {(z + w,zw):\z\<l,\w\<l} (3.2) 

of C 2 . They are relevant to the 2x2 spectral Nevanlinna-Pick problem because, for a 
2x2 matrix A, 

r{A) < 1 o (tr A, det A) G Q 

and 

r(A) < 1 & (tr A, det A) G T. (3.3) 

Accordingly, if F is an analytic 2x2 matrix function on O satisfying r(F(X)) < 1 for all 
A G B then the function (tr F, det F) belongs to the space Hol(B, T) of analytic functions 
from B to r. A converse statement also holds: every cp G Hol(B, T) lifts to an analytic 
2x2 matrix function F on B such that (tr F, detF) = tp and consequently r(F(\)) < 1 
for all A G B (HJ Theorem 1.1]. The 2x2 spectral Nevanlinna-Pick problem can therefore 
be reduced to an interpolation problem in Hol(B, T). There is a slight complication in 
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the case that any of the target matrices are scalar multiples of the identity matrix; for 
simplicity we shall exclude this case in the present paper. 
The relation (|3.3p scales in an obvious way: for p > 0, 



r(A) < p (tr A, det A) G p ■ T 

where 

P- {s,p) d = (ps,p 2 p) and p ■ T d = {p ■ (s,p) : (s,p) G V}. 
The following result is a refinement of [U Theorem 1.1]. 

Proposition 3.1. Let Ai,...,A n be distinct points in B and let W\,...,W n be 2 x 2 

matrices, none of them a scalar multiple of the identity. The following two statements are 
equivalent. 

(1) There exists a rational 2x2 matrix function F, analytic in B, such that 

F(X j ) = W J forj = l,...,n 

and 

supr(F(A)) < 1; (3.4) 

Aeo 

(2) there exists a rational function h G Hol(B, Q) such that 

h(\j) = (trWj.det Wj) for j = 1, . . . , n, (3.5) 
and h(0) is relatively compact in Q . 

Proof. (1)=>(2) Let F be any function with the properties described in (1) and let ro be 
the supremum in the inequality (|3.4p . Let / = (tr,det) o F; then / G Hol(D,<5), / is 
rational and /(B) C ro • T C Q, and so /(B) is relatively compact in Q. 

(2)=^(1) Let h = (hi, ha) be as in (2). Since Q = Uo<t<i * ' ^ an( ^ h(H) is relatively 
compact in Q there exists t G (0, 1) such that /i(B) C t ■ Q. 

Since the 2x2 matrices W\, . . . , W n are not scalar matrices, they are 'nonderogatory', 
that is, their rational canonical forms have only one block, or alternatively, they are similar 
to companion matrices. Hence there exist nonsingular matrices Pi, . . . ,P n such that 



Wj = Pr 



det Wj tr Wj 



Pj, for j = l,...,n. 



Pick a polynomial matrix P(A) such that P(\j) = Pj for each j and P(A) is nonsingular 
for every A G B. The matrix function 



F(A) = P(A)- 1 



1 

-/» 2 (A) hi(X) 



P(A) 



is rational and analytic in B. Its characteristic polynomial is z 2 — h\(X)z + /12(A) and so, 
since hip) Ct-T 

supr(P(A)) <t<l. 

□ 



Certain simple rational functions play a central role in the analysis of T. 
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Definition 3.2. The function <3? is defined for (z,s,p) G C 3 such that zs 7^ 2 by 

• (.,.,,) -5EZi—i. + ^Z^)i. (3.6) 
2 - zs z 1 - 

In particular, $ is defined and analytic onBxT (since |s| < 2 when (s,p) G T), 
extends analytically to (A x T) \{(z, 2z, z 1 ) : z G T}. See [5] for an account of how $ arises 
from operator-theoretic considerations. The 1-parameter family &(uj, •), lo G T, comprises 
the set of magic functions of the domain Q. The notion of magic functions of a domain is 
explained in [8], but for this paper all we shall need is the fact that 

$(D x T) C A 

and a converse statement: if w € C 2 and ^(z,^)! < 1 for all z G B then w E T; see for 
example [3 Theorem 2.1] (the result is also contained in [U Theorem 2.2] in a different 
notation) . 

A T -inner function is the analogue for Hol(D, V) of inner functions in the Schur class. 
A good understanding of rational T-inner functions is likely to play a part in any future 
solution of the finite interpolation problem for Hol(D,r), since such a problem has a 
solution if and only if it has a rational T-inner solution (for example, |14l Theorem 4.2] or 
Theorem 18.11 below) . 

Definition 3.3. A T-inner function is an analytic function h : B — > T such that the radial 
limit 

lim h(r\) e bF (3.7) 

i — >i— 

for almost all AgT, where bT denotes the distinguished boundary ofT. 

By Fatou's Theorem, the radial limit (|3,7p exists for almost all A G T with respect to 
Lebesgue measure. The distinguished boundary &r of Q (or T) is the Silov boundary of 
the algebra of continuous functions on T that are analytic in Q. It is the symmetrisation 
of the 2-torus: 

bT = {(z + w, zw) : \z\ = \w\ = 1}. 
4. An example 

Here is an example of a robust stabilization problem in which all the signals p, q, z, w, y 
and u in Figure 1 are two-dimensional. The theory in Section [2] reduces the robust sta- 
bilizability of this example to a two-point spectral Nevanlinna-Pick problem, which is 
amenable to a precise analysis. 

Example 4.1. Let G = \_Gij\^ be a stabilizable plant such that, for some a G C and 
c^O, 

6163 61 + gb 3 

Lrl1 [ah + fh ^ + 6163]' 

G13 = diag{l,6i/fe 3 }, 

G31 = diag{l, 63/61}, 

61/63" 
63/61 



G, 



33 



(4.1) 



where f, g, 6 7 are as in equations (12.2f) and (|2.3p . 
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Find the values of a, c for which there exists a robustly stabilizing controller for the 
system £ of Figure 1 with respect to the uncertainty set A^o of equation (I2.5p . and describe 
the set of all robustly stabilizing controllers. 

G is unstable, having poles at 1 and 3 in the right halfplane. There do exist stabilizable 
plants G whose corner blocks are as in equations ()4. 1 1) ; examples will be given later. It 
will also transpire that the plant G of Example 12.11 is one such G with a = 10, c = 2. 

Proposition 4.2. There exists a robustly stabilizing controller in Example \A.l\ if and only 
if 

\c\ < 1 —=. (4.2) 

4-2^/3 

Proof. Begin with a doubly coprime factorization of G33. The functions 

"0 b{ 

b 3 

-bi 



N = N 



X 
X 



f 



9 

-63 



9 

-hi 
f 



M = diag{6i,6 3 } 
Y = 



M = diag{&3,M, 



Y 



63 9 

f h 

bi 9 
f b 3 



(4.3) 



belong to RH°° and satisfy equations (|2.10p and (|2.1ip . In the notation of Proposition 
1231 



61 ' 

ab 3 cb^ 



T X = G u + G 13 YMG 31 

T 2 = G 13 M = hi 
T 3 = MG 31 = b 3 L 

Proposition ^. 2l now asserts that £ is robustly stabilizable with respect to the uncertainty 
set Ai if and only if there exists Q € RH°° such that 



X(oo) — N(po)Q(oo) is nonsingular and 

supr((Ti-6i63Q)(s)) < 1- 

sew 

Let k denote the Cayley transform 



(4.4) 
(4.5) 



Ah 



1 + A 
1 - A 



The zeros 1, 3 of the scalar functions T2, T 3 £ H correspond under k 1 to the points 0, \ 
respectively in D. 

Lemma 4.3. Let Ai, A2, C £ ^- There exists a rational function h £ Hol(B, Q) such that 
h(\\) = (C)0)i ^(^2) = ( — C: 0) an d h(P>) is relatively compact in Q (4.6) 
if and only if 

ICI <d(Ai,A 3 ) 

where d denotes the pseudohyperbolic distance on D. 
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Proof. Let 5g,Cg be the Lempert function and Caratheodory distance on Q respectively. 
By definition of the Lempert function, for z\,Z2 G Q there exists h G Hol(B, Q) such that 
h(Xi) = z\ and /i(A 2 ) = ^2 if and only if Sg(zi, z-i) < d(X±, A 2 ). Moreover, if there is such 
an h, then in the case of Q there is a rational h with the same properties. By [7j Corollary 
5.7], <5g = Cg, while (3 Corollary 3.5] gives an explicit formula for Cg. In particular, for 
any (el, 

<5 6 ((C,o),(-c,o)) = c g ((c,o),(-c,o)) = |C|. (4.7) 

=4> Let h G Hol(D,C/) satisfy conditions (|Q]> . Then /i(B) C i • 5 for some f G (0, 1). The 
function 51 = • /i belongs to Hol(B, Q) and satisfies 

g(X 1 ) = (t-\,0), g(X 2 ) = (-t- 1 (,0). 

By definition of Sg, 

^((t- 1 c,o),(-r 1 c,o))<d(Ai,A 2 ). 

Hence, by equation (I4.7j) . 

r i ci<d(Ai,A 2 ), 

and since t < 1 it follows that |£| < (i(Ai,A 2 ). 
<*= Suppose that \(\ < d(Ai, A 2 ). For i G (0, 1], 

d ^ + x t n2 2^|Ax - A 2 | 2 l + t 2 AiA 2 ^ 

— d(t Ai ,tX 2 ) = - — , 2T |2 Re - > 0. 

at |1 — t / A 2 Ai| / 1 — i^AiA 2 

Hence there exists t < 1 such that 

6g(((,0),(-(,0)) = \(\<d(tX 1 ,tX 2 ). 

Consequently there is a rational function g G Hol(B, Q) such that g(tX±) = ((, 0), 5(tA 2 ) = 
(— C,0). Now the function h(X) = g(tX) is rational, is analytic from t~ x B to Q and maps 
Ai,A 2 to (±(,0). Moreover /i(B) C h(B~), a compact subset of Q. □ 

Suppose that G is robustly stabilizable; then we may pick Q G RH°° such that condi- 
tions (|4.4p and (|4.5p hold. Let i 7 = (Ti — T2QT3) o k. This -F is a rational analytic matrix 
function on B such that 



-\a ( v / 3-2)cJ' 



F(o) = r 1 (i) 
F(i) = r 1 (3) 



\ ' 

(2 - v/3)c 



(4.8) 
(4.9) 



and 

supr(F(A)) < 1. (4.10) 

AeO 

The rational function /i = (tr, det) o F belongs to Hol(B, Q), /i(B) is relatively compact in 
Q (see Proposition 13. ID and 

h(Q) = (trF(0),detF(0)) = ((\/3 - 2)c, 0), 

Ml) = (trF(i),detF(i)) = ((2 - V3)c,0). (4.11) 

Hence by Lemma 14.31 

|(V3-2)c| < d(0,|) = |, 
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and therefore 



1 



|c| < = ps 1.866. (4.12) 

4-2^3 

The condition (|4.12p is also sufficient for the existence of a robustly stabilizing controller. 
Suppose it is satisfied. By Lemma 14.31 there exists a rational function h = (/ii,/i2) £ 
Hol(D, C) for some compact subset C of Q such that the interpolation conditions (14.111) 
hold. Let 

/ i (i) = (C,r ? )er. 

Since the matrices Ti(l), Ti(3) in equations (|4.8p are not scalar, there are nonsingular 
matrices Pq,Pi such that PoTi(l)P -1 , FiTi^Pf 1 are companion matrices, that is, 



TU1) 



P 



-i 



1 

(V3-2)c 



P , T 1 (3)=Pf 1 



1 

(2 - V3)c 



Pi. 



Let V be a nonsingular matrix (to be chosen later) and let P be a matrix polynomial such 
that -P(A) is nonsingular for all A 6 B~ and 

P(0) = P , P(^) = Pi and P(1) = V. 

Define P by 



1 

-h,2 h\ 



P. 



Then 



F(0)=Ti(l), F(i) = T 1 (3) and F(l) = V 



-i 



1 

-v C 



v. 



(4.13) 



Now let Q £ RH°° be such that F = (Ti - M3Q) «■ On letting A ->• 1 (and hence 
s = k(A) 00) in this relation we obtain 



V 



-1 



1 

-77 c 



y = F(l) 



(2i - M3QX00) 



and so 

^ I" 1_ 

Q(oo) = 

v ' [a c 
From the equations (14.3|) we have 

Z(oo) = 
iV(oo) = 



F(l) 



1 

a c 



1 

a c 

1 

-77 c 



Q(oo) 



V. 



(4.14) 



9 




-bi 








-63 




f _ 


(00) 






"0 








"0 


1 


&3 





[OO) = 


1 






1 

3 

1 I 



Hence 



X(po) - N(oo)Q(oo) 



_ 1 

3 

-1 

1 

1 



V 



where 



1 

1 

1 

-v C 



1 

a c 
V + Z 



v- 1 



1 

-v C 



V 
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Thus X(oo) — N(oo)Q(oo) is nonsingular provided that 





-V 



VZV 1 is nonsingular. 



(4.15) 



Z is not a scalar matrix, hence it is similar to its companion matrix 

1 

-(f + c+Ia) -2- 

Replacement of Z by its companion form in equation (|4.15p shows that there is a Q € RH C 
such that X (oo) — iV(oo)Q(oo) is nonsingular if V can be found such that 







+ V 







.(io +c+ i a) _ 2 



V" 



(4.16) 



9 1 " 1 3 l 

is nonsingular. A suitable V can be constructed whatever the values of £, 77, o and c, and 
so Q 6 RH°° with the required properties exists. On substituting this Q into the formula 
(|2.12|) for K we obtain the desired robustly stabilizing controller for E. 

□ 



To demonstrate that Example 14.11 is not vacuous we need to show that the plant G is 
stabilizable for a suitable choice of the second block row and column of G. Choose 2x2 
matrix functions R12, R21, R22, R23, R32 £ RH°° arbitrarily and let 

"0 g n 

fh/h 

"0 563/61 




G12 
G21 
G22 
G23 



R12 + 

R21 + -R23 



32 



R22 + R 



R 



23 



G 



32 



63 
61 



23 
6l 

b 3 

1 



gM 
f/h 
1 



R 



32 



R 



32- 



G is now fully specified. If we take a = 10, c = 2 and all the equal to I then we obtain 
the plant G of Example 12. 1L 

To show that G is stabilizable write down a right-coprime factorization of G over RH°° . 
It can be checked that G = hrh~ x where n = r 



*j J ' 



nn 



6163 



«^3 + fh cb^ + 6163 

^21 = R21, fl22 = R22, fl23 = R23 



n\2 = R12, n 13 = bil, 



H31 



1 

/6 3 



"-32 



m 



9 
/ 



32- 



n 33 = iV 



61 
&3 



/ 0' 
0/0 

Ei E 2 M 
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and 



Clearly n, rh G RH C 



where 



and 



E 








o 



9 

f o 



R, 



32- 



Moreover h, rh are right coprime, since 
xrh — yh = I 





' I 





0" 




"0 





0" 


X = 





I 





y = 













£31 


—R32 


X 










Y 



X31 



h 
f 





h 



We claim that the controller K = FX" 1 stabilizes G (and hence G is stabilizable) . Since 
the functions defined in equations (|4.3p satisfy the equation (|2.1ip we have MX — NY = I 
and so X,Y are right coprime. According to \20\ Section 4.2, Theorem 1], K stabilizes G 
if and only if the function 

"0" 



m 



1' 



[0 I]h X 
is invertible in RH°° or equivalently, since [0 / 



m 



I 














/ 














/ 


Y 


C31 


G32 


G33 


X 



[G31 

-1 



G32 G33I rh, if and only if 



(4.17) 



It may be verified that 

















/ 














/ 


Y 


G31 


G32 


G33 


X 



I 






















YCG31 


YCG32 


I + YCG33 


-YC 


—CG31 


—CG32 


—CG33 


C 



where 

C = (X - G33Y) 
Hence K stabilizes G if and only if 

/ 0' 

7 

-M~ l E\ -M' l E 2 M~ l 

I 



diag{6 3 ,6i}. 



10 

1 

YCG31 YCG32 I + YCG33 -YC 

— UU31 — UU32 — OCt33 G 

It is a matter of straightforward calculation to check that this matrix product is 

I 

/ 

63 



e RH C 



f b 3 

-63I 



'^32 

-R32 



X 

-N 



-Y 
M 



/^-SYNTHESIS AS A SEMIDEFINITE PROGRAM 



15 



which belongs to RH°° . Hence G is stabilizable. 

The solution to Example 12.11 now follows: it is the special case in which a = 10, c = 2 
and all the Rij = I. Since Example 14. 1 1 admits a robustly stabilizing controller if and only 
if the inequality (14. 12^ holds, it is not robustly stabilizable when c = 2. 

Proposition 4.4. The answer to the question posed in Example 12.11 is no. There is no 
robustly stabilizing controller for Example 12.11 with uncertainty set Ai o- 

It is striking that in Example 14.11 the set of robustly stabilizing controllers does not 
depend on the five functions Rij. The proof of Proposition 14.21 contains the following. 

Proposition 4.5. The robustly stabilizing controllers K for the system £ of Example \A.l\ 
when \c\ < 1/(4 — 2\/3) are the functions of the form 

K=(Y- MQ){X - NQ)" 1 = (X- QN)~ l (Y - QM) 

where the functions Y,M,X,N,X,N,Y and M are given by equations (14. 3D . 

6i 



Q 



ab 3 06^3 



and F is any rational analytic 2x2 matrix function on D satisfying the conditions 
and (|4.10p and such that 

l 



Fl) + 



r 1 

i-a -1- 



is nonsingular. 



Note that F o k^ 1 is proper for any choice of rational F. 

Because Example 14.11 is constructed so that the resulting model matching problem has 
scalar T 2 and T3 and only two interpolation nodes, the function theory of T is adequate 
for a full analysis of the robust stabilization problem. For the general robust stabilization 
problem, even in the case of the uncertainty set Ai 0, currently known theory of T does 
not suffice to decide the solvability of the corresponding model matching problem. The 
remainder of the paper develops an alternative approach that leads to the criterion in 
terms of a quadratic semidefinite program given in Theorem 1 1 . 1 L 



5. Duality between Hol(D, T) and £2 
The Schur class of the bidisc will be denoted by £2: 

S 2 = f Hol(D 2 ,B-). 

A strategy for the 2x2 spectral Nevanlinna-Pick problem is as follows. 

(1) Reduce to an interpolation problem in Hol(D,r) as in Proposition 13.11 

(2) The magic functions ^(z, •) induce a duality between Hol(D, T) and a subset of £2- 

(3) Use Hilbert space models for 1S2 to obtain a necessary and sufficient condition for 
solvability. 

For the second step, observe that since ^(D x T) C B~, if h = (s,p) G Hol(D,r) then 
the function 

(*, A) ^ h(X)) = 2ZP } X) ~*[ X) for z, A G B 

2 — zs{\) 

belongs to £2- The simple observation that <1> induces a correspondence between Hol(D, T) 
and a subset of £2 underlies the results in this paper. The correspondence was first 
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developed in [3], where a realization theorem for Hol(B, T) was proved and some examples 
were calculated. In this section we answer the question: 

which subset of 52 corresponds to Hol(B,r)? (5-1) 

If h = (s,p) G Hol(B, r) then, for any fixed A £ D, the map 

is a linear fractional self-map f(z) = |fr^ of D with the property "6 = c". To make 
the last phrase precise, say that a linear fractional map / of the complex plane has the 
property "6 = c" if /(0) ^ oo and either / is a constant map or, for some a, 6 and d in C, 

f( z ) = for all z£CU {oo}. 

JK ' bz + d 1 ' 

The following is an easy calculation. 

Proposition 5.1. If f is a non-constant linear fractional transformation then f has the 
property "b = c" if and only if /(0) ^ oo and 

r\z) = - 1 /or alheCU {oo}. 



Here is an answer to the question (15. ip . 

Proposition 5.2. Let G be an analytic function on B 2 . There exists a function h G 
Hol(B, r) such that 

G(z, A) = §(z, h(X)) for all z, A G B (5.3) 

if and only if G £ S2 and, for every A G B, G(-,X) is a linear fractional transformation 
with the property "b = c". 

Proof. Necessity is immediate (note that the relation (|5.3p implies that G(0, A) = — |s(A) 7^ 
00). 

Conversely, suppose that G £ S2 and (?(•, A) is a linear fractional transformation with 
the property "6 = c" for every A G B. Since G(-, A) G ^2 for every A G B, G(-, A) does not 
have a pole at 0, and therefore, by the "6 = c" property, we may write 

a(X)z + b(X) 1 / \ \ (a(A) -6(A) 2 )z , „ ^ _ 
G(z,X) = \> x)z+ \> = 6(A) + V A)z ;; ; for all z,X G B 

for some functions a, 6 on B. Observe that this statement remains true even when G(-, A) 
is constant for some A. Since 

6(A) = G(0,A) 

and 

a(X)z = (b(X)z + 1) (G(z, A) - 6(A)) + b(X) 2 z, 
the functions a and 6 are analytic on B. Let 

s(A) = -26(A), p(X) = a(A) for all A G B, 
and let h = (s,p). Then h is analytic on B and 

g,f 2zp(A) - g(A) a(A)^ + 6(A) 

$(Z ' hW) = 2 - . S (A) = b(X)z + 1 = A) - 
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Since G G 52 we have, for any AeD 

\<&(z,h(X))\ < 1 for all z G O, 

from which it follows that h(X) G T, see [TJ Theorem 2.1 and Corollary 2.2] and [21 
Proposition 3.2]. Thus h G Hol(B,r) has the required properties. □ 



6. The Schur class of the bidisc 

Every function in S 2 has a Hilbert space model [1] . That is to say, if ip G S 2 then there 
exist a separable Hilbert space Ai, a Hermitian projection P on M and an analytic map 
u : B 2 ->■ At such that 

1 - <p(fi)<p(\) = ((I - fj,*pX P )u(X),u(fi)) for all A,/j G B 2 , (6.1) 

where, for A = (Ai, A2) G B 2 , Ap denotes Ai-P + A2(/ — P). This statement is contained in 
the proof of [TJ Theorem 1.12] - see in particular equation (3.11). The triple (M,P,u) is 
called a model of (p. 

The function §>(z, h(\)) has the property that it is linear fractional in z for every A G B. 
A consequence of this property is that the projection P (corresponding to the variable z 
in the defining equation (|6.ip ) has rank one for some model of the function, as we now 
show. 

Denote by 5 2x2 the 2x2 Schur class of the disc, that is, the set of analytic 2x2 matrix 
functions F on B such that ||-F(A)|| < 1 for all A G B. The following result (except for the 
uniqueness statement) is essentially [6l Theorem 1.3]. 

Proposition 6.1. Let h G Hol(B,r). There exists a unique function F = [i^jj G 5 2x2 
such that 

(tr F, det F) = h (6.2) 

and 

Fn = F22, I-F12I = I F21 1 a.e. on T, F\2 is either or outer and -Fi2(0) > 0. (6.3) 
Moreover, for all fj,, A G B and all w, z G C such that 

1 - F 22 {^)w / and 1 - F 22 (X)z / 0, 

F satisfies the identity 



l-$(w,h(fi))$(z,h(\)) = (1 - wz)rf(n, w)7(A, z) 

+ r ) (»,wy(I-F( fJ ,)*F(\))r ] (\,z), (6.4) 



where 



7 (A,z) = (l-F 22 (A)z)^ 1 F 12 (A) and 
z^(X,z) 

rj(X,z) 



1 



(6.5) 



Proof. Let h = (s,p). Consider first the case that s 2 = Ap: then the function F = 
diagj^s, tjs} has the required properties (|6.2p and (|6.3p . and moreover it is the only func- 
tion with these properties. 

In the case that s 2 7^ 4p the H°° function js 2 —pis nonzero and so it has an inner-outer 
factorization, expressible in the form 
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where ip is inner, e G is outer and e (0) > 0. Let 



(6.6) 



Then tr F = s and 
Clearly 



det F = \s 2 - pe G = \s 2 - (\s 2 - p) = p. 



\F l2 \ =e Re 2 G = |F 2 i| a.e. on T, 



F\2 is outer and ^12(0) > 0. Thus F has the properties ()6.2p and (|6.3p . and again it is 
easy to see that F is the only function with these properties. 

We must show that H-FHoo < 1 on D. Let f\ = F12, f 2 = ^21- At almost every point of 

T 



I -F*F 



I 

1 

4 





h 




\\s 


h 


h 


k 




h 


k 



Ap\ 



-2s h - 2f 2 s 



-2f lS - 2f 2 s 



4- s 



Ap\ 



(6.7) 



The diagonal entries of the matrix on the right hand side are non-negative, for if s(A) 
z± + z 2 , p(\) = z\z 2 where z\,z 2 £ B~, then 



4-| S | 2 -| S 2 



Ap\ 



= 4 - \z\ + z 2 \ 2 - \zi - z 2 \ 2 
= A-2\ Zl \ 2 -2\z 2 \ 2 
0. 



(6.8) 



Furthermore, for almost all A € T, 

(-2sf 1 -2f 2 s)(-2f 1 s-2f 2 s) = \2sf 1 + 2f 2 s\ 2 

= 4(| S | 2 |/i| 2 + s 2 hf 2 + s 2 hh + \s\ 2 \h?) 

= 4| S | 2 (|/ 1 | 2 + |/ 2 | 2 ) + 8Re( S - 2 / 1 / 2 ) 

= 2|s| 2 |s 2 -4p| +2Re(s 2 (s 2 -Ap)). 

Hence, for almost all A € T, 

16det(/-F*F) = {A-\s\ 2 -\s 2 -Ap\) 2 -2\s\ 2 \s 2 -Ap\-2Re(s 2 {Ap- s 2 )) 
= 16 + |s| 4 + |s 2 - Ap\ 2 - 8|s| 2 - 8|s 2 - Ap\ + 2|s| 2 |s 2 - Ap\ 

-2|s| 2 |s 2 - Ap\ - 2Re(s 2 (s 2 - Ap)) 
= 16 - 8|s| 2 - 8|s 2 - Ap\ 

+ \s\ 4 + |s 2 - Ap\ 2 - 2Re(s 2 (s 2 - Ap)). 

Note that 



(6.9) 



(6.10) 



|s| 4 + \sf - Ap\ z -2Re{r(s l - Ap)) = 
Thus, for almost all A G T, 

16det(I - F*F) = 16 + 16|p| 2 



Ap)\ 2 = 16|p| 2 . 



8|s 2 -Ap\. 
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Since (s,p) maps D into T, by continuity (s(X),p(X)) can be written as [z\ + ^2,^1^2) for 
some z\,z 2 G B~ , and 



(6.11) 



16det(J-F*F) = 16 + 16|ziz 2 | 2 - 8\zi + z 2 \ 2 - 8\zi - z 2 \ 2 
= 16(1 + \ Zl z 2 \ 2 - \ Zl \ 2 - \z 2 \ 2 ) 
= 16(l-|zi| 2 )(l-|z 2 | 2 )>0. 

The inequalities (|6.8j) and (16. lip show that 

I-F(X)*F(X) > 

for almost all A £ T. Thus F £ S 2x2 . 

We now prove the identity (|6,4p . For A £ D and for z € C such that 1 — ^s(X)z 7^ 0, 



$(z,h(X)) 



2zp(X) - s(X) 
2 - zs(X) 

(p(A) - | S (A) 2 )z 



■i S (A) + 



1 - hs(\)z 



(6.12) 



Next apply a standard type of identity for linear fractional transformations; see [31 
Lemma 1.7]. Let H, U and Y be Hilbert spaces. For any operator 



P 



P21 P22 



h ®u -»• h e y, 



(6.13) 



denote by T-p the linear fractional transformation 

F P (X) = P 22 + P 21 X{I - P u Xr 1 P 12 

defined for any operator X on H such that / — P\\X is invertible. Fp(X), when defined, 
is an operator from U to Y. The following identity of standard type may be verified by 
straightforward expansion. Let P = [Pij] 2 .j = i, Q = [Qij]j-j=i be operators from H ®U 
to H © Y. For any pair of operators X, Y on H such that I — P\\X and / — Q\\Y are 
invertible, 



I-F Q (Y)*F P (X) 



H2 



(i - Y*Q* 11 y 1 (i - y*x) {i - Puxy 1 p 



12 



+ [Q\ 2 (I - Y*Q* u y 1 Y* I] (I - Q*P) 



X(I- PnX)- 1 P 12 
I 



(6.14) 



In the light of the definition (|6.6p of F, for A £ D and for any z £ C such that l—^s(X)z ^ 0, 

1 



i S (A) + / 1 (A)/ 2 (A)z 
±*(A)4 



1 - fs(A)2 

( l s(A) 2_ p(A))z 



and so, by equation (|6.12p . 



1 - ±s(A)z 
F F{x) (z) = -$(z,h(X)). 



(6.15) 
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By the identity (|6.14p . for all /j,AcD and for any i»,zeC such that 1 — ^s(/j,)w / and 
1 — ^s(X)z 0, the expansion 



+ [MiT)(l-w 1 1 s( f ,))- 1 w l)(I-F(»yF(\)) 



= (1 - wz)j(h,w)j{\,z) 

n . „I-F(jjl)*F(\) „ , 
+(l-/2A)7?(M,w)* v{\ z ) 

holds, where 7 and n are defined by equations (|6.5p . 



z (l - i s (A)z)-Vi(A)' 



(6.16) 
□ 



7. An alternative proof of the realization of $>(z,h(\)) on the bidisc 

In this section we give an alternative proof of Proposition 16. II in the special case that h 
is T-inner. 

Proposition 7.1. Let h = (s,p) : ID — > T be a T -inner function. There exist a Hilbert 
space Ai, an analytic function F : D — > C(C 2 ,Ai) and an outer function 7 G H°° such 
that 

| 7 (A)| 2 = 1 - |s(A)| 2 /4 a. e. on T, 

and for all /j,, A € B and /or any iu,zeC suc/i i/iai 1 — ^s(fj,)w 7^ and 1 — |s(A)z 7^ 0, 
the identity 



1 - $(w, h(jj))${z, h(X)) = (1 - iz) 



+(1 - flX) ( F(A) 



7(A) 



t(m) 



1 — i-zs(A) ' 1 — hws(/j,) 



( 


1 > 








27(A) 




f 1 ) 

wy(p) 1 


V 









(7.1) 



holds. 

Proof. Since is T-inner, by [2J Proposition 3.2(3)], s(X)/p(X) = s(X) and |p(A)| = 1 

for almost all A € T. Thus, for almost all A £ T and for z € C such that 1 — ^s(X)z 7^ 0, 



unw z-\s{X)/p{X) 
${z,h(X)) = — 3 — — — p{X) 



1 - \zs(X) 



1 - \zs{X) 



p(X). 



(7.2) 
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By the identity (|7.2j) . since \p{\)\ = 1 for almost all A G T, for to, z G C such that 
1 — 7}s((i)w 7^ and 1 — ^s(X)z 7^ 0, the expansion 



1 - $(u7,/i(A))$(z,/i(A)) 



w - ^s(X) 



p(A) 



1 - ^s(A) 1 - 5*s(A) 



p(A) 



to 



i a(A ) *-f s (A) 



1 - lws(X) 1 - ^zs(X) 
(l-|,(A)| 2 /4) (1-wz) 



(7.3) 



(1 - i fl (A)w) (1 - i S (A)z) 

holds for almost all A € T. By a theorem of F. Riesz [TO], there exists an outer function 
7 G such that 

| 7 (A)| 2 = 1 - |s(A)| 2 /4 a. e. on T. 
Then, for w, z G C such that 1 — ^s([i)w 7^ and 1 — hs(X)z 7^ 0, 



1 - $(w,h(\))$(z,h(\)) = 

Therefore, for all such w,z G C, 
27(A) 107(A) 



7(A) 



7(A) 



1 - |-zs(A) 1 - |ws(A) 



(1 — toz) a. e. on T. (7.4) 



1 + 



1 - lzs{Xy 1 - \ws(X) 



$(w,h(\))$(z,h(\)) + 
7(A) 7(A) \ 



1 - i yZs(A)' 1 - ^WS(X) 



a. e. on T. 



(7.5) 



Thus the relation (|7.5|) can be expressed by the statement that, for almost all A G T and 
for all z G C such that 1 — ^s(X)z / 0, the Gramian of the vectors 



G C 



(7.6) 



/ 1 

is equal to the Gramian of the vectors 

'$(z,h(\)y 

7(A) 
1— 2 ZS W 

Consequently there exists an isometric operator L\ on the vector space C 2 which maps 
the vectors in equation (j7.6[) to the corresponding vectors in equation (|7.7p . 
Define *(A),AgD, by 



(7.7) 



*(A) 



def 



i S (A) 



7 (A) \s{X) 



(7.8) 
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Clearly ^ is analytic on D. It is easy to check that, for all A £ B and for any z G C such 
that 1 - \s(\)z ^ 0, 

/ 1 



¥(A) 



1 

*7(A) 
1— 2^s(A) t 



2 S ( A ) — ^(xy 
7 (A) i s (A) 

7(A) 



(7.9) 



1— 2«s(A) 

Since the boundary values La of ^ are isometries for almost all A G T, the function ^ is 
in the Schur class £2x2- 

For all A G B and for any w,z £ C such that 1 — is(/i)w; 7^ and 1 — |s(A)z 7^ 0, 

/ 1 




Vl-5«(A) 



7(A) 



1 - h(n))$(z, h{\)) - (1 - tDz) 



7(A) 



mwMf 

7(M) 



(7.10) 



1 — |zs(A) ' 1 — ^ius(/x) / 

Since is in the Schur class S2X2, there exists a Hilbert space M and an analytic 
: B -)• £(C 2 , A4) such that, for all /J, A G B, 

/-W(A) = (l-jiAWf(A). 



The dimension of Ai is equal to the rank of 



J-^(/^)*^(A) 
l-/iA 



Therefore, for all fx, A G B and 



for any to,z£C such that 1 — is(//)u> / and 1 — hs(X)z 7^ 0, 



= (1 - wz) ( 




7(A) 




7(A*) 




- \zs(\y 1- 


- iws(/i) 




( 






1 \ 


- *(M)**(A)) 




1 ) 




»7(c) 












= (1 - ute) ^ 




7(A) 




7(m) 




- ±2S(A)' 


1 - 


" 3^s(/^) 



+(1 - /2A) ( F(A) 



/ 


1 > 








27(A) 




1 1 1 


Vi 


~2 ZS ( X )/ 







(7.11) 

□ 
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Remark 7.2. Let h = (s,p) : B — > T be a T-inner function. The relation between the 
analytic function defined by (|7.8p from Proposition 17.11 and the analytic function F 
defined by (16. 6p from Proposition 16.11 is the following. Recall that, for A 6 D, 

p(A)-^ 2 (A)' 



¥(A) 



clef 



i*(A) 



7(A) 

7 (A) |s(A) 
where 7 is an outer function in such that 

| 7 (A)| 2 = l-|s(A)| 2 /4 

Since (s,p) is T-inner, s(A) = s(A)p(A) for almost all A G T, by |2J Proposition 3.2(3)]. 
Therefore s 2 (X)p(X) = |s(A)| 2 for almost all A G T. Hence 



a. e. on T. 



l7(A)f 



1 



| S (A)| 2 /4) 2 = |p(A)-, 2 (A)/4| 2 



a. e. on T. 



Thus, for 



h 



p 



is 2 







and f 2 = 7, 



the functions /1, / 2 G -fP°, /1/2 = s 2 /4 - p and 

|/i(A)| = |/ 2 (A)| a. e. on T. 

Therefore 





h 




"-1 0" 


h 


k 




1 



2 7 



-1 
1 



7 2 s 

8. Criteria for the solvability of spectral Nevanlinna-Pick problems 

The following result (in combination with Proposition 13. ip contains Theorem II .![ 

Theorem 8.1. Let n > 1, let Ai,...,A n be distinct points in B and let (sj,pj) G T 
for j = 1, . . . , n. Let zi,z 2 , z% be distinct points in B. The following five conditions are 
equivalent. 

(1) There exists an analytic function h : B — > T satisfying 

= ( s j,Pj) f° r J = 1, • • • , n; (8.1) 

(2) there exists a rational T-inner function h satisfying equations (JHTTJ) ; 

(3) there exist positive "in-square matrices N = [A^jj.]™^ e k=1 of rank at most 1 and 
M = [Migjk]™f = i £ ^=1 such that, for 1 < i, j < n and 1 < £, k < 3, 



1 



2z £ pi 



ZiSi 



Z k Sj 



(1 - z e z k )Ni£j k + (1 - \i\j)Muj k ; 



1.2) 



(4) there exist positive in-square matrices N = [Ni£j k ]™f = 1 ik=\ °f ran k at most 1 and 



tl7. 



such that 



( 2z e pi - s 
\ 2- Z£Si 



^^kPj Sj 

2 - z k Sj 



> 



[(1 - Z£Z k )Ni£ tjk ] + [(1 - XiX^Mujk) ; 



$.3) 
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(5) the semidefinite program 

min (triV) 2 - tr(iV 2 ) 
subject to the linear matrix inequality (|8.3p and the positivity conditions 



M=[M l£ , jk ]lf =1Ak=l >0, 
is feasible, attains its minimum and has value zero. 
Proof. (2)=>(1) and (3)=^(4) are trivial. 

(3)=>(2) Suppose N, M as described exist. Since N is positive and of rank 1 there exist 
scalars ^j k for j = 1, . . . , n, k = 1, 2, 3 such that 

Likewise, since M > there exist a Hilbert space .M of dimension at most 3n and vectors 
Vj k G such that 

Thus the relation f|8.2j) can be expressed by the statement that the Gramian of the vectors 

-^(z k ,Sj,Pj)\ 

7jk \€C 2 (BM, j = l,...,n, k = 1,2,3, (8.4) 
v jk J 
is equal to the Gramian of the vectors 

Zk-Tjk] eC 2 ®M, j = l,...,n, fc = 1,2,3. (8.5) 

Consequently there exists a unitary operator L on the finite-dimensional vector space 
C 2 @ M which maps the vectors in the expression (|8.5|) to the corresponding vectors in 
the expression (|8.4|) . Write L as a block operator matrix 



L 



A B 
C D 



where A, D act on C 2 , Ai respectively. Then for j = 1, . . . , n, fe = 1,2,3 we obtain the 
pair of equations 

f-$(zk,8j,pj)\ =A ( 1 \ BXjVj k 
Vjk = C {zXk) +DXjVjk - 

From the second of these equations, 



v jk 

and hence 



(8.7) 
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Let 



*(A) = A + BX(I - DX)~ l C 



a(A) 6(A) 
c(A) d(A) 

Since L is unitary and A4 is finite-dimensional, ^ is a rational 2x2 inner function, and 
hence 

h d = (trtf.dettf) (8.9) 

is a rational T-inner function. 

We claim that h satisfies the interpolation conditions (I8.ip . By equation (|8.7p . 

-*(ak,*j,Pi)\ _ \ / 1 \ = ( a(Aj) + 6(A j )z fc 7 jfc N 
Tjfc / jJ \Zk7jk) \c(Xj) + d(\ j )z k j jk/ 

for j = 1, . . . ,n, k = 1, 2, 3. Eliminate 7^-^ from these two equations to obtain 

<&{z k , Sj,pj) = -a(Xj) - b{\j)z k {\ - d(X j )z k )~ 1 c(X j ). 

That is to say that, for each j G {1, . . . , n}, the linear fractional maps 

agree at three distinct values of z € D. It follows that the two maps are identical, which 
is to say that 

a(Xj) = \sj, b(\j)c(Xj) = \sj -pj, d(Xj) = \sj. 

Hence 

tr^(Aj) = a(Xj) + d{Xj) = Sj 

and 

det *(Aj-) = (ad - bc)(Xj) = \s] - {\s) - Pj ) = Pj . 
Thus h(Xj) = (sj,pj) for j = 1, . . . , n as required. Thus (3)=>(2). 

(4)=>(1) The proof of this statement is similar to that of (3) =>-(2). The only difference is 
that the relation (|8.3[) can be expressed by the statement that the Gramian of the vectors 

J jk eC 2 ©M, j = l,...,n, fc = 1,2,3, (8.10) 

v jk ) 

is less than or equal to the Gramian of the vectors 

Zkljkj eC 2 (BM, j = l,...,n, fc = 1,2,3. (8.11) 

Consequently there exists a contractive operator L on the finite-dimensional vector space 
C 2 © M which maps the vectors in the expression (|8.11j) to the corresponding vectors in 
the expression (|8.10p . Since L is contractive, * defined by (|8.8j) is in the 2x2 Schur 
class, and hence 

h = (tr \[/, det G Hol(D, T). (8.12) 
The proof that h(Xj) = (sj, P j) is unchanged. 
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(1)=>(3) Suppose there exists an analytic function h = (s,p) : D — > V satisfying equa- 
tions (I8.ip , By Proposition l6.lt there exists an analytic function F 



c 



2x2 



such that ll-^lloo < 1 on D and, for all fi, A G O and for any w,z £ C such that 1 — ^s(/j,)w ^ 
and 1 - \s{\)z / 0, 



1- $>(w,h(iJ,))$(z,h(\)) = (l-wz)-f(fi,w)~f(\,z) 

+(1 - fi\)rj(fi, w)* x _ i A V(\ *) 



$.13) 



where 



7 (X,z) = (l-i S (A)«)-Vi(A) and 

By assumption, for the given A-,- £ D, j = 1, . . . , n, 

h(Xj) = (sj,pj) for j = l,...,n. 
Let fj, = Aj and A = Aj, 1 < i, j < n, in f|6. 16j) . For all ii),z£D, 



$.14) 
$.15) 



1 -*(iu, supi^fasj^j) = l-<f>(w,h(\i))<f>(z,h(\j)) 

= (1 - wz)i(\i,w)7(\j,z) 

+(1 - A^AjMAi, to)* J ~ f (A f f (Aj V ^, *)• 



.16) 



1 — AjAj 

Let w = Z£, z = Zk, 1 < £, k < 3 in f|8. 16j) . Since no (sj,pj) is equal to (2z k , zfy for any k 



1 - <b(zt, Si,pi)®(z k , Sj,pj) = (1 - zez k )-f(Xi,ze)^(Xj, z k ) 
+(1 - \i\j)rj(\i, zi) 



^-r){Xj,z k ). 



1 — AjAj 



.17) 



Since Ili^Hno < 1 on D, the matricial kernel 



1-/2A 

is positive. Hence the positive 3n-square matrices 

■in, 3 def 



of rank at most 1 and 



r){Xi,z t ) 



J-FiXiTFjXj 
1 — AjAj 



n,3 

i,j=l,e,k=l 



«,3 



i,j=X,i,k=l 



satisfy, for 1 < i, j < n and 1 < £, k < 3, 



$.18) 



1 - $(zt,Si,pi)$(z k ,Sj,pj) = (1 - ztz k )Nu,jk + (1 - KXj)M iejk . 
Therefore (1)=K3). 

For the equivalence of (4) and (5) we require a simple observation involving exterior 
powers of matrices. Recall that for any r x s matrix A = [<Hj\ 1 the second exterior power 
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/\ 2 A is an x (2) matrix whose entries are the 2x2 minors of A. It follows that A has 
rank at most 1 if and only if /\ 2 A = 0. 

Lemma 8.2. (1) For any self-adjoint matrix A, 

2trf\ 2 A = (trA) 2 -tr{A 2 ). 
(2) If A is a positive matrix then 

rank,4<l (tr A) 2 - tr(A 2 ) = ^ (tr A) 2 - tr(A 2 ) < 0. 

Proof (1) Let A = [a id ] . Since A = A* , 

tr(A 2 ) = tr(A*A) = £ la,,) 2 = £4 + 2^ k/, 

ij i i<j 

and so 

2tr /\ 2 -4 = 2^(aiiOjj - |ajj| 2 ) 

= 2 ^ OiiO^- - I tr(,4 2 ) - 4 ) 
= (Y,a ll ) 2 -tr(A 2 ) 

i 

= (trA) 2 -tr(A 2 ). 

(2) Let ^ > 0. Then also /\ 2 A > (if A = 5*5 then by the Cauchy-Binet formula 
A 2 ^ = A 2 (B*B) = (/\ 2 B*)(/\ 2 B) = (A 2 BT(/\ 2 B) > 0). Thus 

rank.4 < 1 f\ 2 A = 

^tr/\ 2 yl < 

4^ (tr.4) 2 -tr(,4 2 ) < 

^ (trA) 2 -tr(A 2 ) = 0. 

□ 

(4)«=>(5) The statement in (4) that there is a feasible pair (N, M) for the LMI ([Q]) 
with N of rank 1 means, in view of Lemma 18.21 that there is a feasible pair for which 
(triV) 2 - tr(iV 2 ) = 0. Since (triV) 2 - tr(7V 2 ) > whenever N > 0, it follows that the 
program in (5) is feasible and has value 0. The argument is reversible, and so (4)^=>(5). 

Thus statements (1) to (5) are all equivalent. 

□ 

Remark 8.3. (1) A natural choice of the points zi,Z2,z% is —1,0,1. This choice is not 
permitted in the theorem as stated above, which requires the to belong to D. However, 
the same proof works for z^ G B~ provided that the denominators 2 — z^Sj in equation 
(18. 2p are nonzero, which is so provided that no (sj,pj) is equal to (2zfc, zi) for any k. 
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(2) The matrix [l — <&(z£, Sj,pj)$(%, Sj,pj)] in statements (3) and (4) of Theorem 18.11 is 
positive if and only if the points (sj,pj) are all equal and lie on the variety s 2 = 4p. 

To see (2) study the proof of (3)=>(2) with 7^ = 0,va = 0. 

(3) The objective function (triV) 2 — tr(iV 2 ), though non-negative and quadratic, is not 



concave. For example, for N 
on (^-dimensional space, 



XI > on n-dimensional space, since /\ I is the identity 



(tr(A/)) 2 -tr ((A/) 2 ) = 2tr (A 2 /) 



A 



which is not concave on The objective cannot be expected to attain its minimum at 

an extreme point of the feasible region. 

The next theorem relates the criterion of Theorem 18.11 to the //-synthesis problem. It 
also includes, for comparison, a statement of another criterion obtained previously in [6]. 



. . , W n be 2x2 complex 



Theorem 8.4. Let Ai, . . . , A n be distinct points in B and let Wi, 

matrices, none of them a scalar multiple of the identity. Let Sj = tr Wj, pj = det Wj for 
each j and let z\,zi,z% be any three distinct points in B. The following four conditions 
are equivalent. 

(1) There exists an analytic 2x2 matrix function F in 



such that 



F&j) = Wj for j = l, 



,n 



and 



5.19) 
$.20) 



r(F(A)) < 1 for all A G ID; 

(2) there exists an analytic function h : B — > T such that 

h(\j) = (trWj.det Wj), j = 1,2, ...,n; (8.21) 

(3) there exists a bounded analytic 2x2 matrix function F inO such that conditions 
(I8.19|) and (18.20j) are satisfied, and in addition, both eigenvalues of F(X) have 
modulus 1 for all X £ T; 

(4) there exist 61, ... , b n , c±, . . . ,c n £ C such that 



I - 






2 s j 


1-, 


\,A 


i 



> 



i.22) 



and 



1, 



,n. 



.23) 



Proof. By Theorem 1.1 and Main Theorem 0.1], conditions (1), (2) and (4) are equiva- 
lent. The equivalence (2)4^(3) follows from (1)^(2) of Theorem 18. II and the fact that, for 
a 2 x 2 matrix A, both eigenvalues of A have modulus 1 if and only if (tvA, det ^4) S bT. 
□ 



Remark 8.5. It is obvious that, under the hypothesis that none of the Wj is a scalar 
multiple of the identity matrix, the conditions (l)-(4) of Theorem 18.41 and the conditions 
(l)-(4) of Theorem 18.11 are equivalent. 
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Remark 8.6. In [6] Main Theorem 0.1], the authors proved that (1) and (4) of Theorem 
18.11 are equivalent under a genericity condition: none of the Wj is a scalar multiple of the 
identity matrix. In [1 H Theorem 1.1] H. Bercovici removed this genericity condition and 
replaced condition (4) by a very similar one. 



9. Matricial formulations of the solvability criterion 



There are more matricial ways of expressing the solvability criteria of Theorem 18.11 
Here are some of them. Just for this section we shall denote by I n the n x n identity 
matrix. 

Theorem 9.1. Let Ai, . . . , A n be distinct points in B and let W\, . . . , W n be 2x2 complex 
matrices, none of them a scalar multiple of the identity. Let Sj = tr Wj, Pj = det Wj for 
each j. Let z\, Z2 and 2:3 be distinct points o/B~ such that no (sj,pj) is equal to (2zk, zf.) 
for any k. 

Let 3n-square matrices X, Z and A be defined by 



X 



1 



/ 2z e pj - Sj \ 2z k pj 



\ 2 — z e si 



Z^Sj 



n,3 



i,j=l,£,k=l 



A = diagjAi}™: 3 , 



1 n,3 



=1) 
= 1' 



Z = diag{*tf} i=lj _ 

The following conditions are equivalent. 

(1) There exists an analytic 2x2 matrix function F in 



such that 



F(Xj) = Wj for 



1. 



, n 



and 



r(F(A)) < 1 for all A e D; 
(2) there exist positive 3n-square matrices N, M such that rank N < 1 and 

X > N — Z*NZ + M- A* MA; 



(9.1) 

(9.2) 
(9.3) 



(9.4) 
(9.5) 
(9.6) 



(3) the same as (2) but for the replacement of > by =; 

(4) there exist a positive 3n-square matrix M , a 1 x 3n vector 7 and a matrix P of 
type 3n x 2 such that 






X 



> 



h 
P 





hn 





7* 




1 

Z*7* 



M 



7 

1 Z 

- A* MA 



h 




p* 



(9.7) 



(5) the same as (4) but for the replacement of > by 

(6) the semidefinite quadratic program 

(trAO 2 - tr(A^ 2 ) 



mm 



subject to the conditions N > 0, M > and the linear matrix inequality (|9,6p is 
feasible and has value 0. 

Note that in N, M, A and Z the rows are indexed by the pair (i, I) and the columns by 
the pair (j, k), where i and j run from 1 to n, and I and k run from 1 to 3. 



30 



JIM AGLER, ZINAIDA A. LYKOVA AND N. J. YOUNG 



Proof. The equivalences (1)<^(2)<£4>(3)<^>(6) are just reformulations of (1)<£4>(3)<^(4)44>(5) 
of Theorem 18. 1[ 

(2)=>(4) Suppose (2). Since rankA^ < 1 and N > there exists a 1 x 3n vector 7 such 
that N = 7*7. Consider the Schur complement identity 



" A 


B 




B* 


D 





h 

B*A~ 1 I- 



3?( 



A 

D — B*A- l B 



h A- l B 

hn 



(9.8) 



where A, D are of types 2 x 2, 3n x 3n respectively. Choose 



,4 



-1 
1 



B 



7 

1 Z 



D = M — A* MA. 



The identity (|9,8p becomes 



-1 





7 













1 












* 

7 


Z*7* M 


- A* MA 










1 


" 






1 










1 










1 







-7* Z*7* 


hn 










M- A* MA + 7*7 



(9.9) 





1 





-7 







1 


1 z 













Let 

Thus equation (|9.9p is 

-1 
1 

7* Z*7* 



p = -P*^" 1 = [7* -z*j*] e c 



3nx2 



/ 2 



7 

7^ 
M - A* MA 

-1 
1 



(9.10) 








M - A* MA + 7*7 - Z*j*'yZ 



h 




p* 



On pre- and post-multiplying by the inverses of the first and third matrices on the right 
hand side and using the relation (|9.6jl we obtain the relation (19. 7h 



h 

P hn 



-10 7 
1 7Z 

7* Z*7* M - A* MA 



J 2 P* 

hn 



-1 







"-1 


0" 


1 





< 


1 








M - A* MA + 7*7 - Z* 7 * 7 Z 







X 



Hence (4) holds. Thus (2)=>(4); the proof that (3)=>(5) is almost identical. 
(4) =^(2) Suppose (4). The inequality (19. 7p can be written 

> 



h 

~P hn 



A 
X 



h 




p* 

hn 



A B 
B* M - A* MA 
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It follows that 







< 



A B 
B* M - A* MA 



Hence P* 



A -AP* 
-PA PAP*+X^ 

-AP* - B 

t -PA - B* PAP* + X - M + A* MA 

-AB and 

< PAP* + X-M + A* MA 

= B*A 3 B + X-M + A* MA 



[7* Z*7*] 



"-1 0" 




7 


1 







+ X - M + A* MA 



= -7*7 + Z* 7 * 7 Z + X - M + A* MA 
and so (2) holds with A" = 7*7. Again, (5)=>(3) is proved in much the same way. 



□ 



Relaxation of the condition rank N < 1 in (2) of Theorem 19. II yields a necessary condi- 
tion for the solvability of a spectral Nevanlinna-Pick problem in the form of the feasibility 
of an LMI. The following statement is immediate from Theorem 19.11 



Corollary 9.2. In the notation of Theorem 19.11 if the spectral Nevanlinna-Pick problem 
(|9.4p - (|9.5p is solvable then there exist positive 3n-square matrices N and M such that the 
inequality (|9.6p holds. 

In fact the existence of positive N and M such that the inequality (|9.6I) holds is equiv- 
alent to the existence of 92 £ £2 such that 

ip(z£,\j) = $(zi,Sj,pj), for t = 1,2,3, j = 1,. .. ,n. 

Since <p{-, A) need not be linear fractional, we cannot derive an h G Hol(D,r) from ip. 

10. Construction of all interpolating functions 

Theorem 18.11 gives us a criterion for the solvability of the interpolation problem 

find h G Hol(D,r) such that h(Xj) = (sj,pj) for j = 1, . . . , n. (10-1) 

The proof of the theorem contains a description of a process for the derivation of a solution 
of the problem (fTOTTj) from a feasible pair (N, M) for the LMI (pUj) or jOD with rank N < 
1. The process can be summarized as follows. 

Procedure SW 

Let zi,Z2, Z3 and Xj, Sj,pj be as in Theorem 18. 11 Let N, M be positive 3n-square matrices 
such that rank N < 1 and the LMI $9M holds. 



■ 1 n,3 

i,j=l,£,k=r 



(1) Choose scalars 7^ such that A^ = [7^7^ 

(2) Choose a Hilbert space M. and vectors £ such that M = [(fjfc, 

(3) Choose a contraction 



1 n >3 
M\i,j=l/,k=V 



A B 
C D 



C 2 (BM^C 2 ®M 
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such that 



for j = 
(4) Let 

for A G 





h(X) = (tr, det)(A + BX(I - DX)~ l C) 



(10.2) 



(10.3) 



Then h G Hol(ID), T) and h(Xj) = (sj,pj) for j = 1, . . . , n. 

The purpose of this section is to show that this procedure in principle yields the general 
solution of the problem (|10.ip , provided that one can find the general feasible pair (N, M) 
for the relevant LMI with rankiV < 1. 

Proposition 10.1. Every solution of a T -interpolation problem arises by Procedure SW 
from a solution (N, M) of the corresponding LMI with rank N < 1 . 

Proof. Let Zk,Sj,pj be as Theorem 18.11 and let h G Hol(B, T) satisfy h(Xj) = (sj,pj) for 
j = 1, . . . , n. We must produce a pair of positive matrices (N, M) that satisfy the LMI 
(|9.6p such that Procedure SW, when applied to (N, M) with appropriate choices, produces 
h. 

By Proposition 16.11 there is a unique F G 5 2x2 such that h = (trF,detF), F 1 
F22, I-F12I = I-F21 1 a.e. on T, F12 is either or outer and -Fi2(0) > 0. Then F\\ = F22 = 
and Proposition 16.11 asserts further that if 



2*' 



7(*,A) 
n{z,X) 
J 



1 



F 12 (A) 

-± s (xy 
1 

1 



and 



I 
then 

1 - <S>(w,h(fi)Mz, h(X)) = (1 - 7Dz)^K7I) 7 (z, A) + v (w, n)*J (I - F(ti)*F(X)) Jn(z, A) 
for all z, w,X,fJ,€ D. Since F G 5 2x2 , the map 



(A,Ai) ^ J 



J-F(/i)*F(A) 
1 - /iA 



J 



is a positive 2x2 kernel on D, and so there is a Hilbert space % and an analytic map 
U : D £(C 2 ,H) such that 



J- 



gg(A) 
1 -/2A 



J = U{p)*U{X) 



for all A, G B. Then 



1 - <&K /i(/x))$(z, /»(A)) = (1 - twMw, M)7(^, A) + (1 - AA)r?(«;, fi)*U(/j,)*U(X)r](z, A). 

(10.4) 
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In particular, when w = zi, \i = Aj, z = Zk, A = Xj, 



1 - <&(zi, h(si,Pi)®(z k , Sj,Pj) = (1 - Z£Zkh(Z£, \i)j(kZ, Xj) 

+ (1 - XiXj) (U(Xj)r](z k , Xj),U(Xi)r](zi, Xi)) n 
for i,j = l,...,n,£,k = l,2, 3. Thus the 3n-square matrices 

M = [(U(Xj) V (z k , X,), U(Xi) V ( Zi , Xi)) H ] 

satisfy the LMI (|9.6p (and even the matrix equation (|8.2p ). and rankiV < 1. We may 
therefore apply Procedure SW to (N,M). In steps (1) and (2) choose 

7jk = j{zk,Xj), M = H, vjk = U(Xj)r](zk,Xj). 
By virtue of the relation (|l(J.4p the Gramian of the vectors 

1 \ 

zj(z,X) gC 2 ©%, z,XeO, (10.5) 
XU(X) V (z,X)J 

is equal to the Gramian of the vectors 

/-*(*, /i(A))\ 

j(z, A) g C 2 e7i, z,AeD. 
\U(X)ri(z,X)J 

Hence there exists an isometry Lq on the subspace of C 2 © T~L spanned by the vectors (110,51) 
such that 

/ 1 \ /-$(z,h(X))\ 
L \ zj{z,X) = 7 (z,A) (10.6) 
\XU(X)r,(z,X)J \U(X) V (z,X)J 



g c(c 2 ®n) 



for all z,AeD. Let 

_ \A B 
\C D 

be any contractive extension of Lq. On specialising equation fj 10.6f> to zt and Xj one 
obtains the relation (|10.2p in step 3 of Procedure SW. One may therefore use L in step 4, 
and so obtain a function h G Hol(D,r) that satisfies h(Xj) = (sj,pj). 
We claim that h = h. By equation (|10.6p . 

U(X)r)(z, A) = C A) ) + DXU(X) V (z, A) 

and so, by elimination of i](z, A), 

(~*& ^ A) )=(^ + BA(/--DA)-»C)>,(z,A) 
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for all z,AeD. Now eliminate j(z, A) to obtain 



for all z,AsD. Since 



2 

it follows that 

*n(A) = ±s(A) = * 22 (A) 

and 

*i 2 (A)* 2 i(A) = i s (A) 2 -p(A). 

Hence 

tr^ = s, det^=p 
and therefore h = h as required. □ 

11. Implementation of the solution procedure 

We conclude with some remarks on the practical feasibility of our results for the numer- 
ical solution of a spectral Nevanlinna-Pick problem. Let interpolation points Ai, . . . , A n 
and target matrices Wi, . . . , W n of type 2 x 2 be given. If any of the Wj are scalar matrices 
then the corresponding interpolation conditions can be removed by the standard process 
of Schur reduction, and so we may suppose that all the Wj are nonscalar. Alternatively, 
if some Wj are scalar, one may still reduce to an interpolation problem for Hol(D,r), but 
with interpolation conditions on derivatives |24| ; one could then try to prove analogs of the 
present results for this wider class of interpolation problems (this should not be difficult). 

Supposing, then, that the Wj are nonscalar, let Sj = tr Wj, pj = det Wj. As recalled in 
Proposition 13.11 the problem reduces to the solution of the interpolation problem (jlO.ip . 
Choose zi,Z2,z$ of modulus at most 1 such that no (sj,pj) is (2zk,zf,) for any k (if 
r(Wj) < 1 for each j then one can make the natural choice of —1,0, 1 for the Zk). 

To determine with the aid of Theorem 11.11 whether the problem (jlO.ip with these data is 
solvable we may test the criterion (3) of the theorem. That is, we must ascertain whether 
there exist positive matrices ./V of rank 1 and M satisfying the LMI fll .3 j) in condition (2). 
Existing software packages can reliably determine whether such an LMI is feasible, but we 
do not know an effective way to test whether there is a feasible pair such that rank N < 1. 
The following refinement of Theorem 18.11 at least shows that a search over a compact set 
of pairs (N, M) suffices. 

Proposition 11.1. Let Xj,Sj,pj and zt be as in J7teorem l8.il The T -interpolation problem 

XjeB^(sj, P j)er, j = i,...,n, (n.i) 

is solvable if and only if there exist positive 3n-square matrices N = [Nnjk\if = i i k=i °f 
rank 1 and M = [M^jfc]"'? =1 e k=1 that satisfy the LMI (|9,6p and 
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and 



1 



(11.3) 



(l-|| Si |)(l-I| Si | 

Proof. Sufficiency is contained in Theorem 18. 11 (3)=^(1). To prove necessity, suppose that 
the interpolation problem is solvable. In the proof of Theorem 18.11 (1)=>(3) it was shown 
that the LMI ([S3]) (or equivalently, ([9"U[) ) holds when 

2 r I n,3 

iV = [N.u,jk\i t j = x,t,k=i = t( a *' ^)7(Aj,2 fc ) 



i,j=l,i,k=l 



of rank 1 and 



ji,3 



M = [Mujk)% 



: J-F(A,)*F(A, 
1 — AjAj 



-r)(\j,Zk) 



/i,3 



i,j=i,e,k=i 



(11.4) 
(11.5) 



where each F(Xj) is a contraction, 

j(Xj,z k ) = (1 - 5S i z fc )" 1 /i(Aj), 

"f(\j,Z k )z k 



T](Xj,Z k ) 



1 



and the function /i is in the Schur class. Thus, for j = 1, 

|7(A,-,«fc)| < V , , 
1 2 1 s i I 



, n and k = 1,2, 3, 



(11.6) 
(11.7) 



from which the estimate (|11.3|) follows. Moreover 



NA;,^)lP 



l{Xj,z k )z k 
1 



^ 1 + 



(l-|| Si | 



and therefore 



from which the bound (|1 1.2|) follows 



11^,^)1111^,^)11 



1 — AjAjl 



i-F(KyF(\ 3 )\\, 



□ 



One approach to the finding of a suitable pair (JV, M) would be to use the alternative 
formulation (2) in Theorem 18.11 to minimize the quadratic function f(N,M) = (tr N) 2 — 
tr(iV 2 ) over the feasible region R. By compactness / attains its minimum on the set of 
positive pairs (N, M) that satisfy the LMI (|9,6p and the bounds (|11.3p and (|11.2p . provided 
that this set is nonempty. Proposition 1 1 1 . 11 asserts that the r-interpolation problem ([11. ip 
is solvable if and only if this minimum is zero. 

Since / is positive homogeneous of degree 2 its local minima over R all lie on the 
topological boundary of R, and the gradient of / is linear in N. However, R is a subset of 
real Euclidean space of 18n 2 variables and the boundary of R has a complicated structure. 

Once a feasible pair (JV, M) for the LMI (|9.6p with rank JV < 1 is found, it is a matter 
of straightforward linear algebra to apply Procedure SW in order to calculate a solution h 
of the interpolation problem ([ll.ip . It is routine to find 3n scalars jj k such that Naj k = 
luljk for all i,£,j,k. Likewise, by Cholesky factorization, one can find 3n vectors Vj k in 
some Hilbert space M. such that Mnj k = {vj k ,vn). Because the LMI (|9.6p holds, there 

that satisfies the relation (|10.2j) ; one may then 



is an isometric operator matrix 



C D 
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define an interpolating T-inner function h by equation (|10.3p . In principle it is simple 
linear algebra to find A, B, C and D when M is chosen to be finite-dimensional. 

As Proposition 110.11 shows, the above procedure when applied to the general feasible 
pair (N,M) for the LMI (19. 6ft with rankiV < 1 yields all possible interpolating functions. 
For numerical implementation one would naturally take M finite-dimensional, and then 
the resulting function h will be rational. A slight modification of Proposition 110.11 shows 
that all rational interpolating functions are obtainable by this procedure. 

An important question is whether the results of this paper furnish an improvement 
(for the special case of the spectral Nevanlinna-Pick problem) on existing l D-K iteration' 
methods, as for instance in the Matlab mu-analysis toolbox [23] . We leave this question 
for future exploration. 
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